试论我国硫铁矿山企业存在的问题及对策   总被引：1，自引：1，他引：0  

黄国庆 《化工矿产地质》2004,26(3):188-192

我国硫铁矿山企业存在着制约发展的国际市场竞争剧烈、债务负担不合理、缺乏公平竞争环境与条件等不利的外部因素,且存在历史遗留的一些内部问题。今后企业发展需在提高经营者素质的同时,还应在争取国家政策扶持、积极推进企业改革、实行“债转股”、分离企业办社会的职能,进行产品结构和产业结构调整等方面制定切实的具体措施加以解决。  相似文献   

一类非线性奇异边值问题解的存在性     

白锦东 《中国海洋大学学报(自然科学版)》1998,(4)

证明一类非线性奇异边值问题解的存在性定理,改进了Ｄｕｎｎｉｎｇｅｒ和Ｋｕｒｔｚ的结果,构造具体的方程说明存在性定理中的条件不能去掉,这表明［４］中一结果不真。  相似文献   

Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data          下载免费PDF全文

Ji Cai  Bülent Tezkan  Yuguo Li 《Geophysical Prospecting》2018,66(8):1602-1624

Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.  相似文献   

Estimation of a non-stationary prior covariance from seismic data     

Rasmus Bødker Madsen  Thomas Mejer Hansen  Henning Omre 《Geophysical Prospecting》2020,68(2):393-410

Non-stationarity in statistical properties of the subsurface is often ignored. In a classical linear Bayesian inversion setting of seismic data, the prior distribution of physical parameters is often assumed to be stationary. Here we propose a new method of handling non-stationarity in the variance of physical parameters in seismic data. We propose to infer the model variance prior to inversion using maximum likelihood estimators in a sliding window approach. A traditional, and a localized shrinkage estimator is defined for inferring the prior model variance. The estimators are assessed in a synthetic base case with heterogeneous variance of the acoustic impedance in a zero-offset seismic cross section. Subsequently, this data is inverted for acoustic impedance using a non-stationary model set up with the inferred variances. Results indicate that prediction as well as posterior resolution is greatly improved using the non-stationary model compared with a common prior model with stationary variance. The localized shrinkage predictor is shown to be slightly more robust than the traditional estimator in terms of amplitude differences in the variance of acoustic impedance and size of local neighbourhood. Finally, we apply the methodology to a real data set from the North Sea basin. Inversion results show a more realistic posterior model than using a conventional approach with stationary variance.  相似文献   

Imaging of elastic seismic data by least-squares reverse time migration with weighted L2-norm multiplicative and modified total-variation regularizations     

Zhiming Ren  Zhenchun Li 《Geophysical Prospecting》2020,68(2):411-430

Least-squares reverse time migration has the potential to yield high-quality images of the Earth. Compared with acoustic methods, elastic least-squares reverse time migration can effectively address mode conversion and provide velocity/impendence and density perturbation models. However, elastic least-squares reverse time migration is an ill-posed problem and suffers from a lack of uniqueness; further, its solution is not stable. We develop two new elastic least-squares reverse time migration methods based on weighted L2-norm multiplicative and modified total-variation regularizations. In the proposed methods, the original minimization problem is divided into two subproblems, and the images and auxiliary variables are updated alternatively. The method with modified total-variation regularization solves the two subproblems, a Tikhonov regularization problem and an L2-total-variation regularization problem, via an efficient inversion workflow and the split-Bregman iterative method, respectively. The method with multiplicative regularization updates the images and auxiliary variables by the efficient inversion workflow and nonlinear conjugate gradient methods in a nested fashion. We validate the proposed methods using synthetic and field seismic data. Numerical results demonstrate that the proposed methods with regularization improve the resolution and fidelity of the migration profiles and exhibit superior anti-noise ability compared with the conventional method. Moreover, the modified-total-variation-based method has marginally higher accuracy than the multiplicative-regularization-based method for noisy data. The computational cost of the proposed two methods is approximately the same as that of the conventional least-squares reverse time migration method because no additional forward computation is required in the inversion of auxiliary variables.  相似文献   

Towards large-scale stochastic refraction tomography: a comparison of three evolutionary algorithms     

Keurfon Luu  Mark Noble  Alexandrine Gesret  Philippe Thierry 《Geophysical Prospecting》2020,68(2):536-552

The main goal of this study is to assess the potential of evolutionary algorithms to solve highly non-linear and multi-modal tomography problems (such as first arrival traveltime tomography) and their abilities to estimate reliable uncertainties. Classical tomography methods apply derivative-based optimization algorithms that require the user to determine the value of several parameters (such as regularization level and initial model) prior to the inversion as they strongly affect the final inverted model. In addition, derivative-based methods only perform a local search dependent on the chosen starting model. Global optimization methods based on Markov chain Monte Carlo that thoroughly sample the model parameter space are theoretically insensitive to the initial model but turn out to be computationally expensive. Evolutionary algorithms are population-based global optimization methods and are thus intrinsically parallel, allowing these algorithms to fully handle available computer resources. We apply three evolutionary algorithms to solve a refraction traveltime tomography problem, namely the differential evolution, the competitive particle swarm optimization and the covariance matrix adaptation–evolution strategy. We apply these methodologies on a smoothed version of the Marmousi velocity model and compare their performances in terms of optimization and estimates of uncertainty. By performing scalability and statistical analysis over the results obtained with several runs, we assess the benefits and shortcomings of each algorithm.  相似文献   

First‐arrival traveltime tomography with modified total‐variation regularization          下载免费PDF全文

Wenbin Jiang  Jie Zhang 《Geophysical Prospecting》2017,65(5):1138-1154

First‐arrival traveltime tomography is a robust tool for near‐surface velocity estimation. A common approach to stabilizing the ill‐posed inverse problem is to apply Tikhonov regularization to the inversion. However, the Tikhonov regularization method recovers smooth local structures while blurring the sharp features in the model solution. We present a first‐arrival traveltime tomography method with modified total‐variation regularization to preserve sharp velocity contrasts and improve the accuracy of velocity inversion. To solve the minimization problem of the new traveltime tomography method, we decouple the original optimization problem into the two following subproblems: a standard traveltime tomography problem with the traditional Tikhonov regularization and a L₂ total‐variation problem. We apply the conjugate gradient method and split‐Bregman iterative method to solve these two subproblems, respectively. Our synthetic examples show that the new method produces higher resolution models than the conventional traveltime tomography with Tikhonov regularization, and creates less artefacts than the total variation regularization method for the models with sharp interfaces. For the field data, pre‐stack time migration sections show that the modified total‐variation traveltime tomography produces a near‐surface velocity model, which makes statics corrections more accurate.  相似文献   

Linear inversion via the discrete wavelet transform pseudoinverse          下载免费PDF全文

M.Z. Fernández‐Muñiz  A. Cernea  J.L. Fernández‐Martínez 《Geophysical Prospecting》2017,65(Z1):131-149

Classical least‐squares techniques (Moore–Penrose pseudoinverse) are covariance based and are therefore unsuitable for the solution of very large‐scale linear systems in geophysical inversion due to the need of diagonalisation. In this paper, we present a methodology to perform the geophysical inversion of large‐scale linear systems via the discrete wavelet transform. The methodology consists of compressing the linear system matrix using the interesting properties of covariance‐free orthogonal transformations, to design an approximation of the Moore–Penrose pseudoinverse. We show the application of the discrete wavelet transform pseudoinverse to well‐conditioned and ill‐conditioned linear systems. We applied the methodology to a general‐purpose linear problem where the system matrix has been generated using geostatistical simulation techniques and also to a synthetic 2D gravimetric problem with two different geological set‐ups, in the noise‐free and noisy cases. In both cases, the discrete wavelet transform pseudoinverse can be applied to the original linear system and also to the linear systems of normal equations and minimum norm. The results are compared with those obtained via the Moore–Penrose and the discrete cosine transform pseudoinverses. The discrete wavelet transform and the discrete cosine transform pseudoinverses provide similar results and outperform the Moore–Penrose pseudoinverse, mainly in the presence of noise. In the case of well‐conditioned linear systems, this methodology is more efficient when applied to the least‐squares system and minimum norm system due to their higher condition number that allows for a more efficient compression of the system matrix. Also, in the case of ill‐conditioned systems with very high underdetermined character, the application of the discrete cosine transform to the minimum norm solution provides very good results. Both solutions might differ on their regularity, depending on the wavelet family that is adopted. These methods have a general character and can be applied to solve any linear inverse problem arising in technology, particularly in geophysics, and also to non‐linear inversion by linearisation of the forward operator.  相似文献   

Full waveform inversion using oriented time‐domain imaging method for vertical transverse isotropic media          下载免费PDF全文

Zhen‐dong Zhang  Tariq Alkhalifah 《Geophysical Prospecting》2017,65(Z1):166-180

Full waveform inversion for reflection events is limited by its linearised update requirements given by a process equivalent to migration. Unless the background velocity model is reasonably accurate, the resulting gradient can have an inaccurate update direction leading the inversion to converge what we refer to as local minima of the objective function. In our approach, we consider mild lateral variation in the model and, thus, use a gradient given by the oriented time‐domain imaging method. Specifically, we apply the oriented time‐domain imaging on the data residual to obtain the geometrical features of the velocity perturbation. After updating the model in the time domain, we convert the perturbation from the time domain to depth using the average velocity. Considering density is constant, we can expand the conventional 1D impedance inversion method to two‐dimensional or three‐dimensional velocity inversion within the process of full waveform inversion. This method is not only capable of inverting for velocity, but it is also capable of retrieving anisotropic parameters relying on linearised representations of the reflection response. To eliminate the crosstalk artifacts between different parameters, we utilise what we consider being an optimal parametrisation for this step. To do so, we extend the prestack time‐domain migration image in incident angle dimension to incorporate angular dependence needed by the multiparameter inversion. For simple models, this approach provides an efficient and stable way to do full waveform inversion or modified seismic inversion and makes the anisotropic inversion more practicable. The proposed method still needs kinematically accurate initial models since it only recovers the high‐wavenumber part as conventional full waveform inversion method does. Results on synthetic data of isotropic and anisotropic cases illustrate the benefits and limitations of this method.  相似文献   

火山Mogi模型形变反演的虚拟观测法及方差分量估计          下载免费PDF全文

王乐洋  温贵森 《地球物理学报》2020,63(5):1775-1786

针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘(Total Least Squares,TLS)联合反演的虚拟观测法,并使用方差分量估计(Variance Components Estimation,VCE)方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值.  相似文献